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ABSTRACT 

We construct the equation of state (EOS) of dense matter covering a wide 
range of temperature, proton fraction, and density for the use of core-collapse 
supernova simulations. The study is based on the relativistic mean-field (RMF) 
theory, which can provide an excellent description of nuclear matter and finite 
nuclei. The Thomas-Fermi approximation in combination with assumed nucleon 
distribution functions and a free energy minimization is adopted to describe the 
non-uniform matter, which is composed of a lattice of heavy nuclei. We treat 
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the uniform matter and non-uniform matter consistently using the same RMF 
theory. We present two sets of EOS tables, namely E0S2 and E0S3. E0S2 is 
an update of our earlier work published in 1998 (EOSl), where only the nucleon 
degree of freedom is taken into account. E0S3 includes additional contributions 
from A hyperons. The effect of A hyperons on the EOS is negligible in the low- 
temperature and low-density region, whereas it tends to soften the EOS at high 
density. In comparison with EOSl, E0S2 and E0S3 have an improved design of 
ranges and grids, which covers the temperature range T = 0.1-10^'^ MeV with the 
logarithmic grid spacing A logio(7'/[MeV]) = 0.04 (92 points including T = 0), 
the proton fraction range Yp = 0-0.65 with the linear grid spacing AYp = 0.01 
(66 points), and the density range pb = lO^'^-lO^^gcm"'^ with the logarithmic 
grid spacing A log]^o(P-B/[s cm~^]) = 0.1 (110 points). 

Subject headings: equation of state — stars: neutron — supernovae: general 
Online-only material: color figures, machine-readable tables 



1. Introduction 

The equation of state (EOS) of dense matter plays an important role in various as- 
trophysical phen omena such as supernova explosions and the formation of neutro n stars 



and black holes ( iJanka et al.l 12007 



Burrows et al.ll2006l : ISumivoshi et al.ll2005L l2009[ ). Sim- 



ulations of core-collapse supernovae cover a wide range of thermodynamic conditions, and 
extremely high density and temperature may be achieved when black holes are formed by 
failed supernovae. The temperature may vary from to more than 100 MeV, the proton 
fraction changes from to around 0.6, and the density can vary from 10^ to more than 
lO^^gcm"^. Clearly, it is very difficult to construct a complete EOS over such a wide 
range of parameters. The information of matter under extreme conditions is far beyond 
our knowledge of nuclear physics from laboratory experiments. Therefore, it is necessary 
to perform a large extrapolation based on a theoretical model that is supported by micro- 



great efforts have been made to study the EOS of nuclear matter {. 


^attimer & Swestv 


1991; 


Schaffner & Mishustin 


1996; 


Weber 


2005; 


Lattimer & Prakash 


200' 


I). However, most of the 



investigations focused on detailed aspects of nuclear matter, which were often restricted to 
the case of zero temperature or high density with uniform distribution of particles. This 
kind of EOS is generally not applicable for use in supernova simulations. So far, there 
exist only two realistic E OSs which are commonly u sed in simulatio ns of core- collapse sup er- 



novae, namely the one by lLattimer fc Swestyl (Il99ll ) and the one by lShen et al.l (Il998bl ). The 
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Lattimer-Swesty EOS is based on a compressible liquid-drop model with a Skyrme force. 
The Shen EOS is based on a relativistic mean-field (RMF) model and uses the Thomas-Fermi 
approximation with assumed nucleon distribution functions in a Wigner-Seitz cell for the 



descrip tion of non-uniform matter. Recently, a Hartree mean-field calculation (jShen et al. 



2010al ) was perform e d for th e Wigner-Seitz cell instead of the Thomas-Fermi approximation 
used by IShen et al.l ( 1l998bl ). The Hartree calculation can incorporate nuclear shell effects, 
but it requires much more computational resources. 



In our earlier work (jShen et al.l ll998al Jbl). we constructed the relativistic EOS for su- 
pernova simulations, which has been widely used in astrophysical simu l ations over the pas t 
decade (jjanka et al.l 120071 : iBurrows et al.l l2006l : ISumiyoshi et al.l l2005l : ISagert et al.l |2009|) . 
The EOS is based on the RMF theory combined with the Thomas-Fermi approximation. 
The RMF theory with nonlinear a and u terms is able to reproduce nuclear matter saturation 
properties and provide a g ood description for both stable and unstable nuclei (ISugahara fc Toki 
1994ai iHirata et al.lll996l ). The Thomas-Fermi approximation in combination with assumed 
nucleon distribution functions and a free energy minimization is adopted to describe the non- 
uniform matter which is modeled as a mixture of a single species of heavy nuclei, alpha parti- 
cles, and free nucleons that exist outside of nuclei, while the leptons can be treated as uniform 
non-interacting particles separately. The RMF results are taken as input in the Thomas- 
Fermi calculation, so the treatments of non-uniform matter and uniform matter in this EOS 
are sufficient to obtain the table in a consistent manner. It would be pr e ferable to treat 



the na ixture of nuclei as recently do ne in iHempel fc Schaffner-BielichI ( I2010l ) , iFurusawa et al. 



( 1201 ll ). and iBlinnikov et al.l (120111 ) for the detailed treatment of electron captures on nuclei 
in supernova core. However, it is beyond the scope of the current update of the Shen EOS 
table, which has been routinely used in astrophysical simulations. 

Toward a more efficient and effecti ve use of the Shen EOS table, we are required to 
make improvements in the EOS given in IShen et al.l (Il998bl ). hereafter referred to as EOSl. 
The main demand is to increase the number of temperature points which is crucial in the 
simulation of core-collapse supernovae. In some cases, the temperature may reach more 
than 100 MeV, so it is encouraged to provide results at some higher temperatures although 
it is probably beyond the applicability of the RMF theory. It is generally believed that at 
sufficiently high temperature and/or density nuclear matter undergoes a phase transition 
to quark-gluon plasma (QGP). Based on experimental data from high-energy heavy ion 
collisions and lattice quantum chromodynamics (QCD) calculations, it is known that the 
critical t emperature for the QCD phase transition is around = 175 MeV for zero baryon 
density ( iGupta et al.l |20 111 ). Therefore, we note that the nuclear EOS at extremely high 
temperature and/or density is not reliable due to the QCD phase transition, although results 
of the RMF theory are provided for practical use in astrophysical simulations. Another 
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suggestion is to use a linear grid for the proton fraction Yp, instead of the logarithmic Yp grid 
used in EOSl, which can add more points in the important region Yp ~ 0.2-0.5 for supernova 
simulations and save memory space by reducing the number of Yp points a t 1^ < 0.1. 



Furthermore, a proton-rich matter may be involved in astrophysical simulations (jPruet et al. 



20051 : iFrohlich et al.l 120061 ). and therefore the information at Yp ^ 0.6 is needed. In EOSl, 



the density grid spacing is only approximately equal at high density, it is now possible to 
perform the calculation with equal grid spacing in the whole range of density. Since it is 
desirable and feasible to make these improvements, we work out a new version of the EOS 
table, hereafter referred to as E0S2, which contains the same degrees of freedom as EOSl. 

In recent years, there has been extensive discu ssion in the l i terature on the influence o f 
non-nucleonic degrees of freedom in dense matter ( iWeberl 120051 : iLattimer &: PrakashI 120071 ) . 
It is generally believed that hyperons appear around twice normal nuclear matter densit y 
in cold neutron star matter (jSchaffner fc Mishustirull996l : IShenI l2002t llshizuka et al.ll2008l ). 
The first hyper on to appear is A that is the lightes t one with an attractive potential in 
nuclear matter ( llshizuka et al.l l2008l : lYue et al.l |2009| ) . S hyperons are now considered to 
appear at a higher density than A, because S hyperons feel a repu lsive potential in nu clear 
matter accor ding to recent d evelop ments in hypernuclear physics ( llshizuka et al.ll2008l ). In 
the work of llshizuka et al.l ( 120081 ). the authors examined the properties of dense matter 
based on an extended RMF model including the full baryon octet. They presented several 
sets of EOS including hyperons for simulations of core-collapse supernovae, which were con- 
nected with EOSl at low density in the simple procedure described in Section 2.3 o f thei r 
paper. The effect of s-wave pion condensation was also examined in llshizuka et al.l ( 120081 ). 
where the authors considere d free thermal pions assuming the pion mass could not be af- 
fected by the interaction. In lOhnishi et al.l (120091 ). the authors examined the possibility of 
s-wave pion condensation in dense matter by using the phenomenological optical potentials 
determined from the pionic atom or pion-nucleus scattering data, and they concluded that 
s-wave pion condensation would hardly take place in neutron stars and especially have no 
chance if hyperons could participate in neutron star matter. The presence of boson con- 
densati on and deconfined quarks in neutron star s has been extensively discussed in many 
works ( ISchaffner fc MishustinI Il996l : IWeberl l2005l : iLattimer fc PrakashI l2007l : lYang fc Shen 
20081 ). It has been suggested that the quark matter may exist in the core of massive neutron 
stars, and the hadr on-quark phase transition can p r oceed through a mixed phase of hadronic 
and quark matter ( lGlendenning]|l992l : IWeberll2005t ILattimer fc PrakashI l2007l : lYang fc Shen 
20081 ). If deconfined quark matter does exist inside stars, it is likely to be in a color supercon- 
ducting phase, and various co l or superconducting p hases have been i ntensively invest i gated 



in re cent years (iBubal 



2010|) and lSagert et al 



a 



2005 



Paulucci et al.l 1201 ll ). In the works of iNakazato et al.l ( 12008 



hood boioh . the authors constructed the EOS tables for simulations 
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of core-collapse supernovae including the hadron-quark phase transition at high density, 
which were connected with EOSl at low density. The hadron-quark phase transition pro- 
ceeded through a mixed phase obtained by the Gibbs conditions for phase equilibrium, 
where the RMF model was used for the hadronic phase and the bag model was adopted 
for the quark phase. Generally, the introduction of non-nucleonic degrees of freedom leads 
to a softening of the EOS and thereby a corresponding reduction in the maximum mass 
of neutron stars. The recent measurement of the Shapiro delay in the radio pulsar PSR 
J1614-2230 yielded a mass of 1.97 ± 0.04 JPemorest et al.lboioi l. Such a high neutron 
star mass provides an important constraint on the EOS at high density and rules out many 
predictions of non-nucleonic components in neutron star interiors. However, it is currently 
difficult to rule out all possible exotica with the 1.97 M0 observation, some theoretical calcu- 
lations inclu ding hyperons and/or qua r ks could be compatible with the observation of PSR 
J1614-2230 Jhattimer fc Prakash|[2007l . boioi Istone et al.ll2007l : IPaulucci et al.lboill ). 



Among these exotic candidates, A hyperons are the most likely to occur in dense mat- 
ter. In a ddition, much more experimental informatio n is now available for A than other 
hyperons ( Takahashi et al.ll200ll : iNakazawa et al.ll2010l ). From the experimental binding en- 
ergies of single-A h ypernuclei, the po tential depth of A in nuclear matter is estimated to be 



around —30 Me V (IShen et al.l 120061). Several recent observations of double-A hypernuclei 



(see Table 4 of iNakazawa et al.l ( 120101 )) indicate that the effective AA interaction sh ould 
be considerably weaker than that deduced from the earlier measurement (lHiyamall2010l ). In 
theoretical studies of A hypernuclei, the RMF theory with nonlinear a and a; terms is a ble to 
provide a reasonable description of single- and double-A hypernuclei ( IShen et al.ll2006l ). The 
influen ce of A hyperons on neu tron star properties has been investigated within the RMF 
model ( ISugahara fc Tokilll994bl ). We would like to examine the effect of A hyperons on the 
EOS for simulations of core-collapse supernovae, while other hyperons such as E and S are 
ignored due to their relatively high threshold densities and lack of available experimental 
data. We construct the relativistic EOS with the inclusion of A hyperons, hereafter referred 
to as E0S3. For the contribution of A hyperons in E0S3, we assume that the equilibrium 
condition /i„ = /xa is valid. In a supernova explosion, the dynamical timescale is of the order 
of milliseconds, which is long enough to establish equilibrium with respe ct to weak interac- 



tions that change strangeness on the timescale of microseconds or less (jSagert et al.l 12009 



2OIOI ). In fact, A hyperons can have a noticeable effect on the properties of matter mainly 
at high density. It is known that weak equilibrium could be achieved at densities above 



10 



12 



g cm ^ where neutrinos begin to be trapped in the core (IKotake et al.ll2006l : iJanka et al. 



20071 ). Therefore, it is justified to use the equilibrium condition /i„ = /xa to determine the 
A fraction of supernov a matter at high den sity where A hyperons have a noticeable contri- 
bution. In the work of llshizuka et al.l (|2008[ ). the authors assumed weak equilibrium among 
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the full baryon octet when they determined the composition of supernova matter at a fixed 
density, temperature, and charge fraction. They examined the appearance of hyperons dur- 
ing the evolution of core collapse and bounce, and found that the effect of hyperons would 
be small in a sphe rical and adiabatic co llapse of a 15 Mq star by the hydrodynamics without 
neutrino transfer ( Ishizuka et ali [20081 ). The hyperons are expected to appear and play an 
important role in t he thermal evolut i on of protoneutron sta rs and the black hole formation 
from massive stars (jPons et al.lll999l : ISumiyoshi et al.l 120091 ). 



We have two aims in this paper. The first is to improve the EOS table according to 
the requirements of the users. In comparison with the earlier version (EOSl), the following 
improvements are made in E0S2 and E0S3 (see Table 1 for details of the comparison). 



• The number of T points is largely increased, the upper limit of T is extended, and 
equal grid spacing for T is used. 

• A linear Yp grid is adopted instead of the logarithmic Yp grid used in EOSl and the Yp 
upper limit is extended. 

• The upper limit of is extended and equal grid spacing is used in the whole range of 
PBi whereas it is only approximately equal at high density in EOSl. 



The finer grids are favorable for better accuracy in the numerical simulations of core-collapse 
supernovae. The wide coverage of the conditions is necessary for numerical simulations of 
astrophysical phenomena including black hole formation, neutron star mergers, and nucle- 
osynthesis. The second aim of this paper is to provide the EOS table with the inclusion of 
A hyperons for the use in simulations of core-collapse supernovae. The difference between 
E0S2 and E0S3 is that the contribution from A hyperons is included in E0S3 when the A 
fraction is larger than 10~^. 

This paper is arranged as follows. In Section [21 we briefly describe the framework to 
calculate the EOS table. We introduce the RMF model with the inclusion of A hyperons 
and explain how to determine the RMF parameters. To make this paper self-contained, we 
give a brief description of the Thomas-Fermi approximation for the non-uniform matter. 
In Section [31 we discuss our results without and with the inclusion of A hyperons. We 
explore the properties of dense matter by examining the phase diagram, compositions, and 
thermodynamic quantities. Section [H is devoted to a summary. In Appendix A, we give the 
definitions of the physical quantities tabulated in the EOS. E0S2 and E0S3 are presented 
in electronic tables and can also be found on several Web sites. In Appendix B, we describe 
various checks made for the EOS tables. 
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2. Model descriptions 



We constr uct the EOS within the RMF framework (jSerot &: WaleckalllQSGt IShen et al. 



199830,1 20061 ) ■ For uniform matter without heavy nuclei formed, the RMF theory can be 
easily used to calculate the properties of matter. For non-uniform matter where heavy nuclei 
are formed in order to lower the f ree en ergy, we adopt the Thomas-Fermi approximation 
based on the work by lOyamatsul (119931 ). The non- uniform matter can be modeled as a 
mixture of a single species of heavy nuclei, alpha particles, and free nucleons that exist 
outside of nuclei. The results of the RMF model are taken as input in the Thomas-Fermi 
calculation, so the treatments of non-uniform matter and uniform matter in this EOS are 
sufficient to obtain the table in a con sistent manner. It would be prefer a ble to treat the 



mixtu re of nuclei as recent l y don e in iHempel fc Schaffner-BielichI ( 120101 ). iFurusawa et al. 



( 120111 ). and iBlinnikov et al.l (l201l[ ) for the detailed treatment of electron captures on nuclei 
in supernova core. However, it is beyond the scope of the current update of the Shen EOS 
table, which has been routinely used in astrophysical simulations. 



2.1. Relativistic mean-field theory 



We adopt the RMF theory with nonlinear a and u j terms to describe homogeneous 
nuclear matter ( ISugahara &: Tokilll994al : IShen et al.ll2006[ ). We study the two cases without 
and with A hyperons. In the RMF approach, baryons interact through the exchange of 
various effective mesons. The exchanged mesons considered in this work include isoscalar 
scalar a nd vector mesons {a and u) a nd an isovector v e ctor meson [p]. In some published 
studies (ISchaffner fc MishustinI Il996l : IShen et al.l l2006l : lYue et al.l 120091 ) , the two hidden- 
strangeness (ss) scalar and vector mesons (a* and 0) were included in a hyperon-rich system. 
It has been foun d that the attrac tion from a* exchange is almost canceled by the repulsion 
from exchange (jShen et al.ll2006l ). Therefore, we neglect the contribution from the exchange 
of these two hidden-strangeness mesons in this work. 



We start with the Lagrangian of the RMF theory including A hyperons. 



RMF 



h,d^ - Ma - g^a - g^^.u^] i/ja 



1 2 2 1 

-m„a 



1 



g2(y 



4 
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where and V'a denote the nucleon and A hyperon fields, respectively, cr, w^, and p'^^ are o", 
and p meson fields with masses mg., m^, and nip. W^'^ and R""^^ are the antisymmetric 
field tensors for and p"'^, respectively. It is known that the inclusion of nonlinear a terms 
is essential to reproduce the properties of nuclei quantitatively and provide a reasonable 
value for the incompressibility, while the nonlinear u term is added in order to reproduce 
the density dependence of the vector part of the nucleon self-energy obta ined in the rela- 
tivistic Brueckner Hartree-Fock (RBHF) theory (jSugahara fc Tokil Il994al). We adopt the 
parameter set TMl listed in Table IH which was determined in ISugahara fc Tokil ( 1l994al ) 
by fitting some ground-state properties of nuclei including unstable ones. With the TMl 
parameter set, the nuclear matter saturation density is 0.145 fm~^, the binding energy per 
nucleon is 16.3 MeV, the symmetry energy is 36.9 MeV, and the incompressibility is 281 MeV. 
The RMF theory with t he TMl parameter set pr ovides an excellent description of nuclear 
matter and finite nuclei ( ISugahara fc Tokil Il994a[ ). and it is also shown to agr ee satisfacto- 



rily w ith experimental data in studies of nuclei with deformed configurations (IHirata et al. 



19961 ). For the parameter s of A hyperons, we use the experimental mass value Ma = 1115.7 
MeV ( lAmsler et al.ll2008l ). Concerning the coupling constants between mesons and A hyper- 
ons, we take g^jg^ = 2/3 based on the naive quark model and q ^lq„ = 0. 621 determined 



by fitting experimental binding energies of single-A hypernuclei (IShen et al.ll2006l l. which 
produce an attractive potential of A in nuclear matter at saturation density to be around 
—30 MeV. It is known that the inclusion of tensor coupling between u and A is important 
to produce small spin-orbit splittings of single-A hypernuclei, but it does not contribute to 
homogeneous matter. The A hyperon is a charge neutral and isoscalar particle, so that it 
does not couple to the p meson. It is shown that these parar neters can reprod uce well the 
experimental data for both single- and double-A hypernuclei (IShen et al.ll2006l ). 



Starting with the Lagrangian ([T]), we deriv e a set of Euler-Lagran ge equations. We 
employ the RMF approximation as described in ISerot fc Waleckal ( 119861 ). where the meson 
fields are treated as classical fields and the field operators are replaced by their expecta- 
tion values. For homogeneous matter, the non-vanishing expectation values of meson fields 
are a = (a), u = (u)^), and p = (p^^). The equations of motion for the meson fields in 
homogeneous matter have the following form: 



P 



mi mi mi 



g^ 



g^^ ij. o„,.\ , fi'< 



m-: 



3p I J, 



m. 



(2) 

(3) 
(4) 
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The stationary Dirac equations for nucleons and A hyperons are given by 



(5) 
(6) 



where stands for the nucleons (N = p or n). = M + g^a and = + Q^f^ are the 
effective nucleon mass and effective A mass, respectively, s denotes the index of eigenstates, 
while and e\ are the single-particle energies. 

In homogeneous matter, baryons occupy single-particle states with the occupation prob- 
ability // [i = p, n, or A). At zero temperature, // = 1 under the Fermi surface, while // = 
above the Fermi surface. For finite temperature, the occupation probability is given by the 
Fermi-Dirac distribution, 



1 



1 



ft 



l + exp[(£|-/.,)/T] 1. 
1 

l + exp[(-£| + /i,)/r] ~7 



exp 



1 



/T 



(7) 



exp 



+ Mf + uA /T 



where i and i denote the particle and antiparticle, respectively, and are the single- 
particle energies. The relation between the chemical potential Hi and the kinetic part of the 
chemical potential z/j is given by 



yUp = t'p + QujUJ + gpP, 

/in = + guj(^ - gpP, 
/iA = Z^A + gf^uj- 

The chemical potential Hi is related to the baryon number density rii as 



Hi 



1 



(9) 
(10) 

fin 



(12) 



where the quantum number s is replaced by the momentum k when we do the integration in 
the momentum space instead of summing over the eigenstates. We denote ub = rip + n„ + tt-a 
as the total baryon number density and Yp = Hp/ ub as the proton fraction in homogeneous 
matter. The coupled equations are solved at fixed ub and Yp. At zero temperature, A 
hyperons appear only at high density when the equilibrium condition /i„ = /iA could be 
satisfied. For finite temperature, a small number of A hyperons may exist at low density 
with /i„ = /xa, and the A fraction increases rapidly at high density. The the r modynamic 
quantities of homogeneous matter have been derived in lSerot &: Waleckal ( 119861 ) . IShen et al. 
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( ll998al Jbl). and ISheru (|2002[ ). so we simply write the expressions here. The energy density of 
nuclear matter including A hyperons is given by 



1 rdkP^k^+Mrift+ff) 

^^ Jo 



i=p,n,A 



the pressure is given by 



i=p,n,A 



dk k' 



k' 







2 2 

2 3 



^yi 2 



3 ^ ^ „i I 1^2,, 2 I ^ ^ , A I ■'■^2^2 



and the entropy density is given by 

1 1"^ 

J2 - dke [-/fln/f-(l-//=)ln(l-/f) 
-/fln/f-(l-/f)ln(l-/^)]. 



i=p,n,A 



(13) 



(14) 



(15) 



2.2. Thomas— Fermi approximation 

In the low-temperature and low-density region, heavy nuclei may be formed in order 
to lower the free energy. We adopt the Thomas-Fermi approximation in combination with 
assumed nucleon distribution functio ns and a fre e ener gy minimization to describe the non- 



uniform matter based on the work by lOyamatsul (119931 ). In this study, we take into account 
the contribution from A hyperons only when the A fraction is larger than 10~^. In 
the region where the heavy nuclei are formed, is quite small, therefore we neglect the 
contribution from A hyperons in the Thomas-Fermi calculation. The non-uniform matter is 
modeled as a mixture of a single species of heavy nuclei, alpha particles, and free nucleons 
that exist outside of nuclei, while the leptons can be treated as uniform non-interacting 
particles separately. For the system with a fixed proton fraction, the leptons play no role in 
the free energy minimization. Hence we mainly pay attention to the baryon contribution in 
this study. 

We assume that each heavy spherical nucleus is located in the center of a charge-neutral 
cell consisting of a vapor of neutrons, protons, and alpha-particles. The nuclei form a body- 
centered-cubic (BCC) lattice to minimize the Coulomb lattice energy. It is useful to introduce 
the Wigner-Seitz cell to simplify the energy of a unit cell. The Wigner-Seitz cell is a sphere 
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whose volume is the same as the unit cell in the BCC lattice. The lattice constant a is defined 
as the cube root of the cell volume, V^eii = = Nb/ub-, where Nb and hb are the baryon 
number per cell and the average baryon number density, respective ly. We define the b aryon 
mass density as pb = Tn^UB with being the atomic mass unit ( lAmsler et al.ll2008[ ). We 
calculate the Coulomb energy using the Wigner -Seitz approximation and adding an energy 
correction for the BCC lattice ( lOyamatsul 119931 ). This energy correction is negligible unless 
the nuclear size is comparable to the cell size. 



We assume the nucleon distribution function n, 



p or n) in the Wigner-Seitz cell 



as 



n.: r 



n 



out 1 



1 - 



u 



(16) 



< r < 

nT\ Ri<r< Keen, 

where r represents the distance from the center of the nucleus and Rceii is the radius of 
the Wigner-Seitz cell defined by Vceii = 4:7iRl^^j3. The density parameters ra™ and ?t,°"* are 
the densities at r = and r > Ri, respectively. The parameters Ri and ti determine the 
boundary and the relative surface thickness of the nucleus. For the distribution function of 
alpha-particle na{r), which should decrease as r approaches the center of the nucleus, we 



assume 



-n 



out 



1 - 



+ n': 



0<r < R. 



n 



out 



'VI 



Rp ^ r < -Rcell, 



(17) 



which could give ^^(r = 0) = and ^^(r > Rp) = 



out 



Here we use the same parameters 
Rp and tp for both proton and alpha-particle distribution functions in order to avoid the 
presence of too many parameters in the minimization procedure. The parameters i?„ and t„ 
may be somewhat different from Rp and tp due to the additional neutrons forming a neutron 
skin in the surface region. For a system with fixed temperature T, proton fraction Yp, and 
baryon mass density pb, there are eight independent parameters among the ten variables, a, 
n™, n°"*, Rn, tn, n^p, n'^^, Rp, tp, and n°"*. The thermodynamically favorable state is the one 
that minimizes the free energy density with respect to these eight independent parameters. 
In principle, the resulting nucleon distribution in the Wigner-Seitz cell would depend on the 
form of the parameterization. It is also possible to determine the nucleon di stribution by a 



self- consistent Th omas-Fermi method without any form of p arameteriza t ion (|Sil et al.ll2001 



Shen et al.ll2010bl ). We have compared, in Figures 1 and 2 of IShen et al.l (l2010bl ). the results 
obtained by the self-consistent Thomas-Fermi method with those using the parameterization 
of Equation f ll6p . and found that they were in good agreement with each other for the cases 
considered. Therefore, Equation f ll6p is considered to be a reasonable form of the nucleon 
distribution in the Wigner-Seitz cell. 
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In this model, the free energy density contributed from baryons is given by 

/ = FceW / 0? = {EceW — T Scew) / Cl^ , 

where the free energy per cell Fcew can be written as 

-^cell = (-Ebulk + Es + Ec) — TSceW = -Fbulk + Eg + Eq- 

The bulk energy -Ebuik, entropy 5*0011, and free energy Fbuik are calculated by 



bulk 



( Un (r) , Up (r) , n„ (r) ) d^r, 



cell 



cell 



"bulk 



(r) , (r) ) (fr, 



cell 



firin (r) 



rir, [r 



Sr. 



cell 



(19) 

(20) 
(21) 
(22) 

Here e ( n„ (r) , rip (r) , n„ (r) ), s ( n„ (r) , np (r) , n„ (r) ), and / ( n„ (r) , np (r) , (r) ) are 
the local energy density, entropy density, and free energy density at the radius r, where 
the system can be considered as a mixed uniform matter of neutrons, protons, and alpha- 
particles. These local densities are the sum of the contributions from nucleons and alpha 
particles. We use the RMF theory described in Section [271] to calculate the nucleon contribu- 
tion, while the alpha-particles are treated as an ideal Boltzmann gas. In general, the number 
density of alpha-particles is quite small, and therefore the ideal-gas approximation is consid- 
ered to be a reasonable approximation for alpha-particles. We note that other treatments for 
alpha-particles have been develo ped and used in nuclear astrophysics (ISumiyoshi fc Ropke 
20081 : [Horowitz fc SchwenkI 120061 ). The free energy density of alpha-particles in the ideal-gas 
approximation is given by 



/„(n,) = -Tn^ [ln(8nQ/n,) + 1] + (4M - S„) 



(23) 



where is the number density of alpha-particles, and we have used the abb reviation ng = 
[MT/ (27r)]^^^. The alpha-particle binding energy Ba = 28.3 MeV is taken from Lattimer fc Swesty 
( I1991I ). We have to take into account the volume of alpha-particle, otherwise the alpha- 
particle fraction would become a large number at high density, where the alpha-particle 
should actually disappear. When we take into account the volume excluded by alpha- 
particles, the free energy densities of nucleons and alpha-particles are given by 



fN{nn,np) 



;i - u)fN{nn,np), 
[1 - u)fo,{h^), 



(24) 
(25) 



where u = riaVa is the fraction of space o ccupied by alpha-particles with the effective volume 
of alpha-particle = 24fm~^ taken from lLattimer fc Swestyl (Il99ll ). We denote the effective 
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number density of neutrons, protons, or alpha-particles as hi = ni/{l — u) {i = n, p, or a). 
The inclusion of the volume excluded by alpha-particles has negligible effect in the low- 
density region, whereas it is necessary for the calculation at high density. 

As for the surface energy term Eg due to the inhomogeneity of the nucleon distribution, 
we take the simple form 



E,= [ Fo I V(n„(r)+np(r)) \U^r. 

J cell 



(26) 



The parameter Fq = 70 MeV fm^ is determined by performing the Thomas-Fermi calculation 
for finite nuclei so as to reproduce the gross propertie s of nuclear masse s, charge radii, and 
the beta stability line as described in the Appendix of lOyamatsul (Il993l ). 



The Coulomb energy per cell Ec is calculated using the Wigri er-Seitz approximation 
with an added correction term for the BCC lattice (jOyamatsulll993l ) 



Ec = l [ e[np (r) + 2n„ (r) - n^] 0(r)ciV + AEc, 

2 Jcell 



(27) 



where 0(r) represents the electrostatic potential calculated in the Wigner-Seitz approxima- 
tion. He is the electron number density of a uniform electron gas {ue = YpUs), and AEc is 
the correction term for the BCC lattice, which is approximated as 



AEc = Q 



BCC- 



(^n 



(2^ 



Here a is the lattice constant, Cbcc = 0.006562 is taken from lOyamatsul (119931 ). and Znon is 
the non-uniform part of the charge number per cell given by 



Rr. 



Anr dr. 



(29) 



Because of the long-range nature of the Coulomb interaction, Ec is dependent on the lat- 
tice type. This dependence has been extensively discussed in lOyamatsul ( 119931 ). The sys- 
tem prefers the BCC lattice because it gives the lowest Coul omb energy. In general, the 
Cou lomb energy is depend e nt on the temperature as discussed in lBravo &: Garcia-Senzl (119991 ) 
and lPotekhin &: Chabrierl (120101 ). The dependence of T in Equation ( 127|) is contained in the 
pa rticle distribution s. We assume that Equation ( I28p . which was derived at zero temperature 



m 



Oyamatsul ( 119931 ). remains valid at finite temperature. 



Since we assume the lattice of nuclei with the Wigner-Seitz approximation, we do not in- 
clude a contribution from the translational energy of heavy nuclei. This contribution is small 
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in general and does not affect the general behavior of the EOS. It may be noticeable only 
within the limited region of temperature, where nuclei behave as gas wi thout being dissoci- 
ated i nto nucleons. We note that this is different from the treatment by iLattimer &: Swesty 
( ll99ll ). who included this term as a minor contribution. 



3. Results 

In this work, we construct the EOS tables covering a wide range of temperature T, 
proton fraction Yp, and baryon mass density pb for the use of core-collapse supernova sim- 
ulations. We present two sets of EOS tables, namely E0S2 and E0S3. E0S2 takes into 
account only the nucleon degree of freedom, while E0S3 includes additional co ntributions 



from A hyperons. In comparison with the earlier version (EOSl) described in IShen et al. 



(Il998bl ). several improvements are made in E0S2 and E0S3 according to the requirements 
of the users. We largely increase the number of T points, and adopt a linear Yp grid in E0S2 
and E0S3 instead of the logarithmic Yp grid used in EOSl. On the other hand, the numerical 
methods are improved to allow the calculation with equal grid spacing for pB, while it is only 
approximately equal at high density in EOSl. For a detailed comparison between the EOS 
tables discussed in this paper, one can see Table 1 in Section [H In principle, the matter at 
extremely high temperatures and densities is beyond the applicability of the RMF theory, 
but we still include the results in these exotic regions since they are sometimes necessary in 
astrophysical simulations. 



3.1. Properties of matter without hyperons 

In E0S2, we present results of matter without the inclusion of A hyperons. For each 
T, Yp, and pb, we determine the thermodynamically favorable state that has the lowest 
free energy in the present model. We perform the free energy minimization for both non- 
uniform matter and uniform matter. Here the phase of heavy nuclei formed together with free 
nucleons and alpha-particles is referred to as non-uniform matter, while the phase of nucleons 
mixed with alpha-particles without heavy nuclei is referred to as uniform matter. For non- 
uniform matter, the minimization procedure is realized by using the Thomas-Fermi method 
which includes eight independent parameters as described in Section [2^21 For uniform matter, 
we perform the minimization with respect to converting two protons and two neutrons into 
an alpha-particle, in which there is only one independent parameter. By comparing the free 
energies of non-uniform matter and uniform matter, we determine the most favorable state 
and estimate the phase transition between non-uniform matter and uniform matter. 



- 15 - 



We first discuss tlie pliase diagram of nuclear matter at finite temperature. It is known 
that tlie density of tlie pliase transition between uniform matter and non-uniform matter 
depends on both T and Yp. The non-uniform matter phase can exist only in the low- 
temperature region (T < 14 MeV). In Figure [H we show the phase diagram in the p^-T 
plane for Yp = 0.1, 0.3, and 0.5. The shaded region corresponds to the non-uniform matter 
phase where heavy nuclei are formed. The dashed line is the boundary where the alpha- 
particle fraction Xa changes between < 10~^ and X^ > 10~^. It is shown that heavy 
nuclei can exist in the medium-density and low-temperature region. The phase of nuclear 
matter at low density is a homogeneous nucleon gas with a small fraction of alpha-particles, 
the heavy nuclei are formed at some medium densities where the system can lower the free 
energy by forming heavy nuclei, and it becomes uniform matter as the density increases 
beyond ~ lO^^'^gcm"^. It is seen that the starting density of the non-uniform matter 
phase depends on T strongly, while the ending density is nearly independent of T. As the 
temperature increases, the density range of the non-uniform matter phase becomes narrower, 
and it disappears completely for T > 14 MeV. In Figure |2l we show the phase diagram in the 
PB^Yp plane for T = 1, 10, and 100 MeV. The shaded region corresponds to the non-uniform 
matter phase, while the dashed line is the boundary between X^ < 10~^ and X^ > 10"^. It 
is shown that the dependence of the boundary on Yp is relatively weak except at very small 
values of Yp. In the case of T = 10 MeV (middle panel), the non- uniform matter phase 
disappears at lower Yp because it is difficult to form heavy nuclei with smaller values of Yp. 
The alpha-particles still exist at T = 100 MeV, and X^ reaches to be more than 10~^ at 
some medium densities for 1^ > 0.1 as shown in the top panel of Figured W e note that 



the pha se diagram of nuclear matter at zero temperature has been discussed in IShen et al. 
(1l998al jb[). 

In the non-uniform matter, nucleon distributions are determined by minimizing the free 
energy density with respect to the independent parameters in the Thomas-Fermi approxima- 
tion. The heavy nuclei are assumed to form a BCC lattice in order to minimize the Coulomb 
lattice energy. In Figure |3l we show the nucleon distributions along the straight line joining 
the centers of the nearest nuclei in the BCC lattice for the case of T = 1 MeV and Yp = 0.3. 
It is found that the nuclei become heavier and get closer with each other as the density 
increases. In principle, there exist free nucleons and alpha-particles outside heavy nuclei at 
finite temperature, but their densities are too small to be observed in Figure |3l We show in 
FigureHlthe same quantities as in Figure[3]but for the case of T = 10 MeV and Yp = 0.3. We 
note that the starting density of the non-uniform matter phase in this case is 10^^'^ g cm~'^, 
so the heavy nucleus does not exist at pb = lO^^gcm"^ for T = 10 MeV which is different 
from the case of T = 1 MeV as shown in the bottom panel of Figure [31 Comparing with 
the top and middle panels of Figure [SJ it is seen that there are much more free nucleons and 
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alpha-particles outside nuclei, while the nucleons inside nuclei get less. This is because the 
role of entropy becomes more important at higher temperature and the free energy can be 
lowered if the nucleons are freed from nuclei. In Figure we plot the nuclear mass number 
A and charge number Z as a function of the baryon mass density pb for Yp = 0.3 at T = 1 
MeV and T = 10 MeV. It is shown that A and Z have relatively weak dependence on pB at 
lower density and they increase rapidly just before the phase transition at pb ~ lO^^gcm"^. 
For the same ps and Yp, the values of A and Z at T = 10 MeV are much smaller than those 
at T = 1 MeV. This is because more nucleons are freed from nuclei at higher temperature, 
and it is eventually impossible to form heavy nuclei for T > 14 MeV in the present model. 

In Figure |6l we show the fraction of neutrons, protons, alpha-particles, and heavy nuclei 
as a function of the baryon mass density pb for Yp = 0.3 at T = 1, 10, and 30 MeV. At low 
density, the matter is a uniform gas of neutrons and protons with a small fraction of alpha- 
particles. The alpha-particle fraction Xa increases with increasing pB, but the formation of 
heavy nuclei at low temperatures causes a rapid drop of Xa, Xp, and X„, which is due to the 
fact that heavy nuclei use up most of the nucleons in non-uniform matter. When the density 
increases beyond ~ lO^^'^gcm"^, the heavy nuclei dissolve and the matter becomes uniform. 
The alpha-particles may exist up to pb ~ 10^^'^ gcm~^ where the volume excluded by alpha- 
particles plays an important role that it is unfavorable to have alpha-particles in the uniform 
matter at such high density. For T > 14 MeV, the heavy nuclei cannot be formed, but there 
are finite values of Xa, especially at medium densities as shown in the top panel of Figure |6l 
The alpha-particle fraction decreases with increasing temperature, and we find that Xa at 
T = 100 MeV is of the order of 10~^ at some medium densities. The alpha-particle fraction 
may be significantly affected if the alp ha-particle binding energy Ba is dependent on the 
density and temperature as discussed in lRopke et al.l ( l2005l ). For simplicity, we neglect this 
dependence and take Ba = 28.3 MeV in the present model. Since we treat the uniform 
matter and non-uniform matter consistently using the same RMF theory, all the resulting 
thermodynamic quantities are consistent and smooth in the whole range. We will discuss 
the thermodynamic quantities in E0S2 and compare with those in E0S3 in the next section. 



3.2. Properties of matter with the inclusion of A hyperons 

In this section, we discuss results of matter with the inclusion of A hyperons given in 
E0S3. It is found that the contribution from A hyperons is negligible in the low-temperature 
and low-density region. In fact, at zero temperature A hyperons appear only at high density 
when the equilibrium condition pn = Pa could be satisfied. At low temperature, such as 
T = 10 MeV, the A fraction Xa is smaller than 10~^ at densities below normal nuclear 
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matter density. For simplicity, we take into account the contribution from A hyperons only 
when Xa is larger than 10~^. In the non-uniform matter phase, Xa is quite small, therefore 
we neglect the contribution from A hyperons in the Thomas-Fermi calculation. 

In Figure [TJ we show the fraction Xi {i = A, n, p, a, or heavy nuclei A) as a function 
of pb for Yp = 0.3 at T = 1, 10, and 30 MeV. For lower temperatures, such as T = 1 
MeV (bottom panel) and T = 10 MeV (middle panel), there is a significant fraction of A 
hyperons only at high density. It is shown that X\ increases with increasing p^, which 
causes a decrease of X^. We note that Xp, which is equal to Yp in the uniform matter at 
high density, has been fixed to be 0.3 in this figure. For T = 30 MeV (top panel), X\ is 
of the order of 10~^ at low density, and increases rapidly at high density. Comparing with 
the top panel of Figure El it is seen that the inclusion of A hyperons does not affect Xa 
and Xp significantly. We find the contribution from A hyperons increases with increasing 
temperature. At T = 100 MeV, Xa is about 12% at low density for the case of 1^ = 0.3. 
On the other hand, Xa decreases with increasing Yp at fixed ps and T, since there are less 
neutrons at higher Yp to realize the equilibrium condition pn = Pa- 

We now discuss the thermodynamic quantities in E0S3 and compare with those in 
E0S2. Here we mainly discuss their properties at high density where there are noticeable 
differences between the results with and without A hyperons. As for the behavior of th er- 
modynamic quantities at low density, one can refer to our earlier work ( Shen et al. 1998al lbl). 



In Figure [HI we show the free energy per baryon F as a function of the baryon mass density 
Pb with Yp = 0.1 and 0.5 at T = 1, 10, and 100 MeV. The results with A hyperons given in 
E0S3 are shown by solid hues, while those without A hyperons given in E0S2 are displayed 
by dashed lines. It is found that the inclusion of A hyperons can lower the free energy and 
this effect increases with increasing pb- However, this effect decreases with increasing Yp 
as shown in Figure [HI which is due to the fact that Xa decreases with increasing Yp. We 
show in Figure [HI the pressure p as a function of pb with Yp = 0.1 and 0.5 at T = 1, 10, 
and 100 MeV. The pressure is calculated from the derivative of the free energy as described 
in Appendix A. Therefore, the effect of A hyperons on the pressure is similar to the one 
observed in Figure [H It is obvious that the inclusion of A hyperons tends to soften the EOS 
at high density. In Figure [TUj we show the entropy per baryon 5* as a function of with 
Yp = 0.1 and 0.5 at T = 1, 10, and 100 MeV. At T = 1 MeV (bottom panel) it is hard to 
see the difference between the results with and without A hyperons, while at T = 10 MeV 
(middle panel) there are small differences at high density. For the case of T = 100 MeV (top 
panel), the effect of A hyperons can be seen in the whole range of density, since Xa reaches 
~ 16% (~ 9%) at low density for Yp = 0.1 (0.5) a t T = 100 MeV. In general, the inclusion 



of A hyperons tends to increase the entropy. In IShen et al.l (ll998al Jbl) , we have discussed 
the effect of the formation of nuclei on the entropy, which has a strong Yp dependence. The 
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behavior of S ai pB < 10^^ gem ^ in the middle and bottom panels of Figure [TO] is due to 
the formation of heavy nuclei in non-uniform matter. 

For neutron star matter at zero temperature, the inclusion of A hyperons leads to a 
softening of the EOS of neutron star matter and thereby a corresponding reduction in the 
maximum mass of neutron stars. We use E0S3 to calculate the neutron star properties 
and find that the maximum mass of neutron stars is about 1.75 Mq, whereas the value 
for EOSl and EQS2 is 2.1 8 M(t^. According to the recent measurement of PSR J1614-2230 



(1.97 M©; iDemorest et al.l ( l2010l )). E0S3 seems to be too soft due to the inclusion of A 
hyperons, while EOSl and E0S2 are compatible with the observation of PSR J1614-2230. 
This is a common difficulty for EOSs with the inclusion of hyperons. So far, there are large 
uncertainties in the properties of matter at high density, which are crucial for determining 
the maximum mass of neutron stars. 



4. Summary 



In this paper, we have presented two sets of the EOS tables (E0S2 and E0S3) covering 
a wide range of temperature T, proton fraction Yp, and baryon mass density pb for the use of 
core-collapse supernova simulations. The difference between E0S2 and E0S3 is that only the 
nucleon degree of freedom is taken into account in E0S2, while E0S3 includes additional 
co ntributions from A hyperons. In comparison with the earlier version (EOSl) described 



m iShen et al.l ( 1l998bl ). several improvements have been made in E0S2 and E0S3 according 
to the requirements of the users. We have largely increased the number of T points, and 
adopted a linear Yp grid in E0S2 and E0S3 instead of the logarithmic Yp grid used in EOSl. 
In addition, we have performed the calculation with equal grid spacing in the whole range 
of density, while it is only approximately equal at high density in EOSl. We have presented 
a detailed comparison between these EOSs in Table 1. 

We have employed the RMF theory with nonlinear a and u terms, which can provide 
an excellent description of nuclear matter and finite nuclei including unstable ones. The 
Thomas-Fermi approximation in combination with assumed nucleon distribution functions 
and a free energy minimization has been adopted to describe the non-uniform matter that 
is modeled as a mixture of a single species of heavy nuclei, alpha particles, and free nu- 
cleons outside of nuclei. The RMF results have been taken as input in the Thomas-Fermi 
calculation, so the treatments of non-uniform matter and uniform matter in this EOS are 
sufficient to obtain the table in a con sistent manner. It would be prefer a ble to treat the 



mixtu re of nuclei as recent l y don e in iHempel fc Schaffner-BielichI (120101 ). iFurusawa et al. 



(j201ll ). and iBlinnikov et al.l (l201l[ ) for the detailed treatment of electron captures on nuclei 
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in supernova core. However, it is beyond the scope of the current update of the Shen EOS 
table, which has been routinely used in astrophysical simulations. 

We have included the contribution from A hyperons in E0S3, when the A fraction is 
larger than 10~^. It is based on the consideration that the A hyperon is most likely to 
occur in dense matter among all possible non-nucleonic degrees of freedom because it is the 
lightest hyperon with an attractive potential in nuclear matter. The potential depth of A in 
nuclear matter is estimated to be around —30 MeV from the experimental binding energies 
of single-A hypernuclei. Several recent observations of double-A hypernuclei suggest that the 
effective AA interaction should be weakly attractive. The RMF theory with nonlinear o and 
bj terms has been extended to includ e hyperons, and i t can provide a reasonable description 



of single- and double-A hypernuclei (IShen et al.ll2006l ). We have not included contributions 



from other hyperons, such as S and S, due to their relatively high threshold densities and 
lack of available experimental data. In E0S3, we have performed a consistent calculation 
including A hyperons for the entire table, whereas the EOS tables of high density matter 
w ere connected with E OSl at low density in the simple procedure described in Section 2.3 



of llshizuka et al.l (120081 ) . so E0S3 is considered to be constructed in a consistent manner for 
the inclusion of A hyperons. We have examined the effect of A hyperons on the properties 
of dense matter. The contribution from A hyperons is negligible in the low-temperature and 
low-density region, whereas it tends to soften the EOS at high density. 

In principle, the matter at extremely high temperatures and densities is beyond the 
applicability of the RMF theory, but we still include the results in these exotic regions since 
they are sometimes necessary in astrophysical simulations. It will continue to be a challenge 
for nuclear physics to provide realistic EOS, which should be supported by microscopic theory 
and available experimental data, for use in astrophysical studies. 
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A. Physical quantities in the EOS tables 

The EOS tables are available on the Web at 
http:/ /physics, nankai. edu. cn/ grzy/shenhong/EOS/ index, html^ 
http://www. rcnp. osaka-u. ac.jp/'^shen/ 
http:/ / user, numazu-ct. ac.jp /'^ sumi/ eos /index, html 

The table under the name "*.tab" is the main EOS table, while "*.tOO" and "*.ypO" are 
those for T = and 1^ = 0, respectively. E0S2 and E0S3 are presented in electronic tables, 
Tables 3-8, in this paper. We present the main EOS table of E0S2 (E0S3) in Table 3 (Table 
6), which is 149 MB (167 MB) in size. Table 4 (Table 7) is the one of E0S2 (E0S3) at T = 
which is 1.6 MB (1.8 MB), while Table 5 (Table 8) is the one for Yp = which is 2.3 MB 
(2.6 MB). The tables are written in the order of increasing T. In the tables on the Web, the 
values of log^g(T) and T (in MeV) are given at the beginning of each block, and the blocks 
with different T are divided by the string of characters "cccccccccccc" . In Tables 3-8 of this 
paper, \og^Q{T) and T are set in columns 1 and 2. We note that the entry for log]^Q(T) is set 
to —100 in the tables for the case of T = 0. 

For each T, we present the results in the order of increasing Yp and pb- The quantities 
in one line of Tables 3-8 are defined as follows. 



1. Logarithm of temperature: logio(^) [MeV]. 

2. Temperature: T [MeV]. 

3. Logarithm of baryon mass density: logj^olPs) [gcm~^]. 

4. Baryon number density: ub [fm~'^]. 

The baryon number density is related to the baryon mass density as pb = mun R with 
rriu = 931.494 MeV being the atomic mass unit taken from lAmsler et al.l (120081 ). 

5. Proton fraction: Yp. 

The proton fraction Yp of uniform matter is defined by 

~ ~ — ; ; n — ' y-^'-) 



where Up, n„, n\, and Ua are the number density of protons, neutrons, A hyperons, and 
alpha-particles, respectively, ub is the total baryon number density. For non-uniform 
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matter case, Yp is the average proton fraction defined by 



Y.-^, (A2) 



where Np and Nb are the proton and baryon numbers per cell given by 

Np = [ [np{r) + 2n^{r)]d\ (A3) 

7 cell 

Nb = [nnir) + np{r)+An^{r)]d^r. (A4) 



J cell 

Here, np{r) and n„(r) are the proton and neutron distributions given by Equation ( fT6ll . 
and na{r) is the alpha-particle distribution given by Equation ( JTTll . Because the A 
fraction is very small in the non-uniform matter phase, we neglect the contribution 
from A hyperons in this case. 

6. Free energy per baryon: F [MeV]. 

The free energy per baryon is defined relative to the free nucleon mass M = 938 MeV 
in the TMl parameter set as 

f 

F = J--M. (A5) 

Ub 

7. Internal energy per baryon: i^int [MeV]. 

The internal energy per baryon is defined relative to the atomic mass unit m„ = 931.494 
MeV as 

^int = (A6) 

Hb 

8. Entropy per baryon: S [fc^]- 

The entropy per baryon is related to the entropy density via 

S = —. (A7) 

Hb 

9. Mass number of the heavy nucleus: A. 

The mass number of the heavy nucleus is defined by 

^= [rinir) + np{r)]Anr'^dr, (A8) 
Jo 

where is the maximum of Rp and which is considered as the boundary of the 
heavy nucleus. 
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10. Charge number of the heavy nucleus: Z. 

The charge number of the heavy nucleus is defined by 

np (r) Anr'^dr. (A9) 

11. Effective nucleon mass: [MeV]. 

The effective nucleon mass is obtained in the RMF theory for uniform matter. In 
the non-uniform matter phase, the effective nucleon mass is a function of space due 
to inhomogeneity of the nucleon distribution, so it is meaningless to list the effective 
nucleon mass for non-uniform matter. We replace the effective nucleon mass by 
the free nucleon mass M in the non- uniform matter phase. 

12. Free neutron fraction: 

The free neutron fraction is given by 

Xn = «"V°"*)/(^BKell), (AlO) 

where Keii = = AttRI^^JS is the cell volume, V""^^ = Keii - 4:TtR\/3 is the volume 
outside the nucleus, n°"* is the number density of free neutrons outside the nucleus, 
and is the average baryon number density. 

13. Free proton fraction: Xp. 

The free proton fraction is given by 

Xp = (n°"V°"*)/(nBKeii), (All) 
where n°^^ is the number density of free protons outside the nucleus. 

14. Alpha-particle fraction: X^. 

The alpha-particle fraction is defined by 

X„ = mj{nBV,eii), (A12) 
where is the alpha-particle number per cell obtained by 

iV„ = / (r) d^r (A13) 



J cell 

and na{r) is the alpha-particle distribution given by Equation (fT7|l . 

15. Heavy nucleus fraction: Xa- 

The heavy nucleus fraction is defined by 

Xa = A/{nBVcen), (A14) 
where A is the mass number of the heavy nucleus as defined in Equation 
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• 16. Pressure: p [MeVfm 

The pressure is calculated through the thermodynamic relation 

p= [nl{dF/dnB)]^^y^. (A15) 

• 17. Chemical potential of the neutron: /i„ [MeV]. 

The chemical potential of the neutron relative to the free nucleon mass M is calculated 
through the thermodynamic relation 

= [ d{nBF)/dnn+A , (A16) 

where Un+A = (1 — Yp) hb is the sum of the neutron and A densities. 

• 18. Chemical potential of the proton: Hp [McV]. 

The chemical potential of the proton relative to the free nucleon mass M is calculated 
through the thermodynamic relation 

f,p^[d{nBF)/dnp]^^^^^^. (A17) 

• 19. Effective A mass: [MeV]. 

The effective A mass is obtained in the RMF theory for uniform matter. We replace 
the effective A mass by the free A mass M\ — 1115.7 MeV when the A hyper on is not 
taken into account. 

• 20. A fraction: X\. 
The A fraction is given by 

ATa = ua/ub, 

where n\ is the number density of A hyperons in the uniform matter phase 
that X\ = is adopted in the non-uniform matter phase. 

B. Checks on the EOS tables 

We have done the following checks on the EOS tables. 

• 1. Consistency of the fractions: 

Xp + X^ + Xa + X^ + Xa^ 1. (Bl) 

• 2. Consistency of the relation between F, -Eint, and S: 

F = Ei^t -TS + rUu-M. (B2) 
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3. Consistency of the thermodynamic quantities: 



V 

F = - Yp) + ^pYp . 



(B3) 



In general, these consistency relations can be satisfied with ir i a few thousandths. The physical 
constants used in this study are taken from lAmsler et al.l (120081 ). 
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Table 1. Comparison between the EOS tables discussed in this paper 







EOSl 


EOS2 


EOS3 


Constituents 


Uniform matter 


n, p, a 


n, p, a 


n, p, a, A 




Non-uniform matter 


n, p, a, A 


n, p, a, A 


n, p, a, A 


T 


Range 


-1.0 < logio(T) < 2.0 


-1.0 < logio(T) < 2.6 


-1-0 < logio(T) < 2.6 


(MeV) 


Grid spacing 


Alogio(T) -0.1 


Alogio(T) =0.04 


Alogio(T) = 0.04 




Points 


32 (including T = 0) 


92 (including T = 0) 


92 (including T = 0) 




Range 


-2 < logioCip) < -0.25 


0<Yp< 0.65 


0<Yp < 0.65 


Yv 


Grid spacing 


Alogio(yp) = 0.025 


AYp = 0.01 


AYp = 0.01 




Points 


72 (including Yp = 0) 


66 


66 


Pb 


Range 


5.1 < logio (ps) < 15.4 


5.1 < logio (Pis) < 16 


5.1 < logio(Pfl) < 16 


(g cm-3) 


Grid spacing 


Alogio(PB) ^ 0.1 


Alogio(pB) = 0.1 


Alogio(pB) = 0.1 




Points 


104 


110 


110 


Reference 




Shen et al. f iggSbl 


This paper 


This paper 



Note. — EOSl has been described in IShen et al] l|l998bl ). while EOS2 and EOS3 are given in this study. EOSl 
and EOS2 take into account only the nucleon degree of freedom, while E0S3 includes additional contributions from 
A hyperons. The non- uniform matter is modeled as a mixture of free neutrons (n), free protons (p), alpha-particles 
(a), and a single species of heavy nuclei (A). The A fraction in non-uniform matter is quite small, so we neglect 
the contribution from A hyperons in the non-uniform matter phase of EOS3. In addition, we add the results of zero 
temperature, this gives a total of 92 points for T in EOS2 and EOS3. 



Table 2. The parameter set TMl for the RMF Lagrangian 



M rric 


rriui rup 




9p 


32 (fm-l) 


93 


C3 


938.0 511.19777 


783.0 770.0 10.02892 


12.61394 


4.63219 


-7.23247 


0.61833 


71.30747 



Note. — 



The masses are given in MeV. 



Table 3. The main EOS table of E0S2 



logio(T) 


T 


loglo(PB) 


»B 




F 




S 


A 


Z 










Xa 


P 


Jin 




(MeV) 


(MeV) 


(g cm-=*) 


(fm-3) 




(MeV) 


(MeV) 


(fcs) 






(MeV) 










(MeVfm-3) 


(MeV) 


(MeV) 


-1.0 


0.1 


5.1 


7.581B-11 


0.01 


-1.524 


6.408 


14.27 


91.37 


29.40 


938 


0.9689 








0.03107 


7.344E-12 


-1.222 


-21.73 


-1.0 


0.1 


5.2 


9.544B-11 


0.01 


-1.502 


6.408 


14.04 


91.39 


29.32 


938 


0.9688 








0.03116 


9.245E-12 


-1.199 


-21.78 


-1.0 


0.1 


5.3 


1.201B-10 


0.01 


-1.480 


6.408 


13.82 


91.56 


29.33 


938 


0.9687 








0.03121 


1.164B-11 


-1.176 


-21.83 



Note. — It covers the temperature range T = 10 ^-lO"^'^ MeV with the logarithmic grid spacing A log^Q (T/[MeV]) = 0.04 (91 points), the proton fraction range Yp = 0.01-0.65 with the linear 
grid spacing AYp = 0.01 (65 points), and the density range PB = 10^"^-10^^ gcm~^ with the logarithmic grid spacing A log^Q (pB/[g cm~^]) = 0.1 (110 points). 
(This table is available in its entirety in a machine- readable form in the online journal. A portion with less digits is shown here for guidance regarding its form and content.) 
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Table 4. The table of E0S2 at T = 



logio(T) 


T 


loglo(PB) 


»B 




F 




S 


A 


Z 










Xa 


P 


Mn 




(MeV) 


(MeV) 


(g cm-3) 


(fm-3) 




(MeV) 


(MeV) 


(fcs) 






(MeV) 










(MeVfm-3) 


(MeV) 


(MeV) 


-100 





5.1 


7.581B-11 


0.01 


-0.2446 


6.261 





96.11 


28.40 


938 


0.9661 








0.03383 


5.848E-16 


3.571E-05 


-24.46 


-100 





5.2 


9.544B-11 


0.01 


-0.2446 


6.261 





96.11 


28.40 


938 


0.9661 








0.03383 


8.889E-16 


4.156E-05 


-24.46 


-100 





5.3 


1.201B-10 


0.01 


-0.2446 


6.261 





96.12 


28.41 


938 


0.9661 








0.03383 


1.328B-15 


4.770B-05 


-24.46 



Note. — It covers the proton fraction range Yp = 0.01-0.65 with the linear grid spacing AYp = 0.01 (65 points), and the density range pB = lO^'^-lO^® gem ^ with the logarithmic grid spacing 
^^ogiQ{pB/[scm-^\) = 0.1 (110 points). 

(This table is available in its entirety in a machine- read able form in the online journal. A portion with less digits is shown here for guidance regarding its form and content.) 



I 

CO 

I 



Table 5. The table of E0S2 for K, = 



l°glo(T) 
(MeV) 


T 
(MeV) 


loglo(P-B) 


(fm-3) 




F 
(MeV) 


^int 
(MeV) 


S 
(kB) 


A 


Z 


(MeV) 


x„ 




Xoc 


Xa 


p 

(MeVfm~^) 


Mn 
(MeV) 


(MeV) 


-100 





5.1 


7.681E-11 





2.281E-03 


6.608 











938 


1.0 











-1.703E-13 


3.663E-06 


-938 


-100 





5.2 


9.644E-11 





2.240E-03 


6.608 











938 


1.0 











-2.099E-13 


4.142E-06 


-938 


-100 





5.3 


1.201E-10 





1.151E-03 


6.607 











938 


1.0 











-1.323E-13 


4.828E-06 


-938 



Note. — It covers the temperature range T = lO""*^— lO"^'^ MeV with the logarithmic grid spacing A logj^Q{T/[MeV]) — 0.04 (92 points including T — 0) and the density range — ID^'"*^ — 
lO^^gcm"^ with the logarithmic grid spacing A log^Q (pg /[g cm~^]) = 0.1 (110 points). 

(This table is available in its entirety in a machine-readable form in the online journal. A portion with less digits is shown here for guidance regarding its form and content.) 



Table 6. The main EOS table of E0S3, which has the same ranges and grids for T, Yp, and pb as the one of E0S2 

given in Table 3 



logio(T) 


T 


loglo(Ps) 






F 




S 


A 


Z 








Xoc 


Xa 


P 


fin 




MX 


Xa 


(MeV) 


(MeV) 


(g cm-3) 


(fm-3) 




(MeV) 


(MeV) 


(fcs) 






(MeV) 










(MeVfm-3) 


(MeV) 


(MeV) 


(MeV) 




-1.0 


0.1 


6.1 


7.681E-11 


0.01 


-1.624 


6.408 


14.27 


91.37 


29.40 


938 


0.9689 








0.03107 


7.344E-12 


-1.222 


-21.73 


1116.7 





-1.0 


0.1 


6.2 


9.644E-11 


0.01 


-1.602 


6.408 


14.04 


91.39 


29.32 


938 


0.9688 








0.03116 


9.246E-12 


-1.199 


-21.78 


1116.7 





-1.0 


0.1 


6.3 


1.201E-10 


0.01 


-1.480 


6.408 


13.82 


91.56 


29.33 


938 


0.9687 








0.03121 


1.164E-11 


-1.176 


-21.83 


1116.7 






Note. — Comparing with EC)S2. two additional columns arc included, which are the effective A mass and the A fraction . 

(This table is available in its entirety in a machine-readable form in the online journal. A portion with less digits is shown here for guidance regarding its form and content.) 
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Table 7. The table of E0S3 at T = 0, which has the same ranges and grids for Yp and as the one of E0S2 given 

in Table 4 



loglo(T) 


T 


loglo(PB) 






F 


Bint 


S 


A 


Z 










Xa 


P 


Mn 




Ml 




(MeV) 


(MeV) 


(g 


(fm-3) 




(MeV) 


(MeV) 


(fcs) 






(MeV) 










(MeVfm~^) 


(MeV) 


(MeV) 


(MeV) 




-100 





6.1 


7.681B-11 


0.01 


-0.2446 


6.261 





96.11 


28.40 


938 


0.9661 








0.03383 


5.848E-16 


3.571B-05 


-24.46 


1115.7 





-100 





6.2 


9.544B-11 


0.01 


-0.2446 


6.261 





96.11 


28.40 


938 


0.9661 








0.03383 


8.889E-16 


4.156E-05 


-24.46 


111,"').7 





-100 





6.3 


1.201B-10 


0.01 


-0.2446 


6.261 





96.12 


28.41 


938 


0.9661 








0.03383 


1.328B-15 


4.770B-05 


-24.46 


1115.7 






Note. — (This table is available in its entirety in a machine- read able form in the online journal. A portion with less digits is shown here for guidance regarding its form and content.) 
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Table 8. The table of E0S3 for Yp — 0, which has the same ranges and grids for T and pb as the one of E0S2 given 

in Table 5 



logio(T) 


T 


loglo(PB) 






F 


^'int 


S 


A 


z 










Xa 


p 








Ml 


Xa 


(MeV) 


(MeV) 


(g cm-3) 


(fm-3) 




(MeV) 


(MeV) 


(kB) 






(MeV) 










(MeVfm" 




(MeV) 


(MeV) 


(MeV) 




-100 





6.1 


7.681B-11 





2.281E-03 


6.508 











938 


1.0 











-1.703E- 


13 


3.553E-05 


-938 


1115.7 





-100 





6.2 


9.544B-11 





2.240E-03 


6.508 











938 


1.0 











-2.099E- 


13 


4.142E-()5 


-938 


1115.7 





-100 





6.3 


1.201B-10 





1.151E-03 


6.507 











938 


1.0 











-1.323E- 


13 


4.828B-05 


-938 


1115.7 






Note. — (This table is available in its entirety in a machine- read able form in the online journal. A portion with less digits is shown here for guidance regarding its form and content.) 
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Fig. 1. — Phase diagram of nuclear matter at 1^ = 0.1, 0.3, and 0.5 (bottom to top) in the 
Pb^T plane. The shaded region corresponds to the non- uniform matter phase where heavy 
nuclei are formed. The dashed line is the boundary where the alpha-particle fraction 
changes between Xa < 10~^ and Xa > 10~^. (A color version of this figure is available in 
the online journal.) 
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Fig. 2. — Phase diagram of nuclear matter at T = 1, 10, and 100 MeV (bottom to top) in 
the pB^Yp plane. The shaded region corresponds to the non- uniform matter phase where 
heavy nuclei are formed. The dashed line is the boundary where the alpha-particle fraction 
Xa changes between < 10~^ and X^ > 10"'^. The non-uniform matter phase does not 
exist at T = 100 MeV (top). (A color version of this figure is available in the online journal.) 
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Fig. 3. — Density distributions of neutrons (solid) and protons (dashed) along the straight 
line joining the centers of the nearest nuclei in the BCC lattice for the case of T = 1 MeV 
and Yp = 0.3 at pb = 10^^, lO^^'^, and lO^^ '^gcm"^ (bottom to top). (A color version of this 
figure is available in the online journal.) 




Fig. 4. — Density distributions of neutrons (solid), protons (dashed), and alpha-particles 
(dotted) along the straight line joining the centers of the nearest nuclei in the BCC lattice 
for the case of T = 10 MeV and Yp = 0.3 at ps = lO^^-^ and lO^^-^gcm"^ (bottom to top). 
(A color version of this figure is available in the online journal.) 
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Fig. 5. — Nuclear mass number A (top) and charge number Z (bottom) as a function of the 
baryon mass density pb for 1^ = 0.3 at T = 1 MeV (sohd) and T = 10 MeV (dashed). (A 
color version of this figure is available in the online journal.) 
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Fig. 6. — Fraction of neutrons (dash-dotted), protons (dotted), alpha-particles (dashed), 
and heavy nuclei (solid) as a function of the baryon mass density pb for = 0.3 at T = 1, 
10, and 30 MeV (bottom to top). (A color version of this figure is available in the online 
journal.) 



-40 - 



■L'.'JIJ.'.UIi.'.'J.li.'.'J.liJ.'±L!J.'J.L'-!JJ.I.UIiJ.UlJJ.'J.L; I ' 



10-' b- 

10-' 
10-= 



10" 

10"-' b- 



10'" ■LLUJi i . i ilL i JilL i JliUJiiLJJiLLlJiLU i jJ ' l y * 
X. 



10-" 
10-^ 
10-^ 
10"' 
10-" 



T=30 MeV 



.LA 



I I 

' I 

' I 

'■ ' I ■ 



T=10MeV| \ / 

X ,' 



X , 

X ' 



i ' . ' -lli ' . ' -ili ' . ' 1 1 1 1 1 1 1 1 1 i jj 1 1 n 1 1 1 1 1 1 1 1 1 1 1 1 



■ / 



. I . . ■ I . . ■ . I . . . . I 



T=l MeV 



.■I ... >i .... I .... I . 



1 1 M 1 1 I h 



5 6 7 



9 10 11 12 13 14 15 16 



iogio(pJ [g/cm^] 



Fig. 7. — Same as Figure El but with the inclusion of A hyperons. (A color version of this 
figure is available in the online journal.) 
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Fig. 8. — Free energy per baryon F as a function of the baryon mass density ps with 1^ = 0.1 
(blue) and 0.5 (red) at T = 1, 10, and 100 MeV (bottom to top). The results with A hyperons 
given in E0S3 are shown by solid lines, while those without A hyperons given in E0S2 are 
displayed by dashed lines. (A color version of this figure is available in the online journal.) 
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Fig. 9. — Same as Figure IH1 but for the pressure p. (A color version of this figure is available 
in the online journal.) 
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Fig. 10. — Same as Figure IH1 but for the entropy per baryon 5*. (A color version of this 
figure is available in the online journal.) 



